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Abstract. This paper discusses evidence for and properties of disks associated to brown dwa rfs in the star- 
forming region p Oph. We selected nine objects from the ISOCAM survey of Bontemps et al. 2001 that have 
detections in the two mid-infrared bands (6.7 and 14.3 /jm), relatively low extinction and low luminosity. We 
present low-resolution near-infrared spectra in the J, H and K bands, and determine for each source spectral type, 
extinction, effective temperature and luminosity by comparing the spectra to those of field dwarfs and to the most 
recent model stellar atmospheres. The results indicate that eight objects have spectral types M6-M7.5, effective 
temperature of 2600-2700 K, one has a later spectral type (M8.5) and lower temperature (about 2400 K). The 
derived extinctions range between Av~2 and 8 mag. The location of the objects on the HR diagram, in spite of 
the uncertainties of the evolutionary tracks for young objects of substellar mass, indicates that all the objects 
are very young and have masses below about 0.08 Mq. The coolest object in our sample has mass in the range 
8-12 Mj (0.008-0.012 Mq). In all cases, the mid-infrared excess is consistent with the predictions of models of 
disks irradiated by the central object, showing that circumstellar disks are commonly associated to young brown 
dwarfs and planetary-mass objects. Finally, we discuss possible variations of the disk geometry among different 
objects, as well as the possibility of using these data to discriminate between various formation scenarios. 
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1. Introduction 

A large number of objects with sub-stellar mass are now 
known, with masses ranging from the hydrogen burn- 
ing limit that divides stars from brown dwarfs (BDs; 
M*<0.075 Mq) to values comparable to the mass of giant 
planets and below the deuterium burning limit (< 0.013 
Mq). Their discovery in regions of star formation has pro- 
voked an intense debate on the formation mechanism of 
such objects. Do they form, as solar mas s star s do, from 
the collapse of a molecular core (Shu et al. |l987| )? Are they 
stellar embryos, whose further growth is prevented by dy- 
namical ejections from small stellar systems (Reipurth & 
Clarke |200l| ; Bate et al. |2002| )? Or are they "planets", 
i.e., objects that form in gravitationally unstable regions 
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of circumstellar disks (Papaloizou & Terquem 2001; Lin 
et al. 1998Q ? Is there a single formation process for all 
substellar objects? What is the lowest mass for the gravi- 
tational collapse mechanism? 

A crucial contribution to this debate is expected from 
studies of the circumstellar disks (if any) associated with 
sub-stellar objects, since different theories make very dif- 
ferent predictions. Disks are a necessary step in any for- 
mation mechanism that involves accretion from a parental 
core. If BDs form from core collapse, they should be asso- 
ciated to disks similar in properties to those found around 
low mass pre- main-sequence stars (T Tauri stars; TTS). 
A prediction of the stellar embryo theory is that the disks 
should be truncated by the ejection mechanism, so that 
they should be small and short-lived. In the "planetary" 
hypothesis, any circumstellar disk should be even less sub- 
stantial. 

In some young BDs, emission in excess of that due to 
the phot osphere has bee n det ected in the nea r (Oa sa et 
al. |1999| ; Muench et al. |200l| ; Wilking et al. |l999D and 



mid-infrared (Persiet al. 200C 



Comeron et al. 199S , 2000 



Bontemps et al. |200l|) , and has been interpreted, by anal- 
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ogy with TTS, as evidence for circumstellar disks. In an 
earlier study (Natta & Testi 2001; Paper I), we discussed 
the properties of three objects in Chamaeleon I for which 
we could find in the literature ground-based spectroscopy 
and photometry as well as ISOCAM measurements at 6.7 



0.5 



and 14.3 /xm (Comeron et al. p000| Persi et al. gOOCj ). One 
of these objects is a bona-fide BD, while the two others are 
close to the threshold between stars and BDs. We found 
that the excess emission was clearly detectable only in the 
mid-infrared, because the stellar photosphere overwhelms 
the disk emission in the three near-infrared bands. The 
observed SEDs are well described by disk models simi- 
lar to those of TTS, assuming that the heating is due to 
irradiation from the central star. 

This first result provides strong support for the idea 
that BDs form like stars, from the contraction of a molec- 
ular core. Hence, we decided to extend our study of disk 
properties to a larger number of substellar mass objects in 
regions of star formation, possibly down to objects of few 
Jupiter masses. With this in mind, we selected a small but 
well defined sample of nine objects in the p Oph region, 
that were detect ed by ISOCAM at both 6.7 and 14.3 x<m 
(Bontemps et al. [2001 ). We obtained near-infrared spectra 
for all of them (see §2), which we used to derive the basic 
parameters of the central objects, namely effective temper- 
ature, luminosity and mass (§3). Because of the adopted 
selection criteria, all of these objects have excess emission 
in the mid-IR. We model the expected disk emission for 
each object and show the results in §4. We discuss the 
implications of our findings in §5 and present conclusions 
in 56. 



2. Observations and data reduction 

2.1. Selection criteria 

We chose nine BD candidates from the sample of Class II 
obj ects d etected at both 6.7 and 14.3 /xm by Bontemps et 
al. (2001). We selected all objects with visual extinction 
less than ^8.5 mag and lumin osity less than ~0.04 L Q 
according to Bontemps et al. ( [2001 ). The first criterion 
ensures the possibility of obtaining high signal to noise 
specta across the entire near infrared range. The low lu- 
minosity was required to increase the chance of selecting 
objects in the range of masses we are interested in. 

The location of the selected objects in the ISOCAM 
color-magnitude diagram is shown in Figure (filled cir- 
cles). We note that all the objects are close to or below 
the completeness limit of the ISOCAM survey, as expected 
for such low luminosities. In colors, our sample span the 
whole range covered by the Class II objects (essentially 
classical TTS). 

Some of the selected sources were known from previ- 
ous studies to be very low-mass objects. In Table 1 we 
give the ISO source number, the J2000 coordinates, other 
designations and references to previous studies. Finding 
charts are provided in Appendix |a|. We will comment on 
the comparison between the literature source parameters 
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Fig. 1. ISOCA M co lor-magnitude plot (adapted from 
Bontemps et al. 2001 ). The symbols are asteriscs for Class 
I sources, open circles for Class II sources, stars for Class 
III sources. Filled circles show the objects in our sample. 
The dotted line indicates the ISOCAM completeness limit. 



and those derived in this paper in Section |Tl]. The results 
on one of the sources in our sample (#033) have been al- 
ready presented in Testi et al. (2002a, hereafter Paper II); 
they have been re-analyzed and reported again here for an 
easier comparison with the rest of the sample. 

2.2. Near-infrared spectroscopy 

Near-infrared spectra for the objects in our sample were 
acquired in the period July 4-9, 2001 at the Telescopio 
Nazionale Galileo (TNG), using the multi-mode Near- 
Infrared Camera Spectrogr aph ( NICS; Baffa et al. 2001 ). 
The Amici device (Oliva 200C), a prism based, high- 
throughput optical clement unique to NICS, was used as 
disperser, coupled with a 0.5" wide, 4.2' long slit; the re- 
sulting effective resolution is approximately A A/ A ^100, 
approximately constant across the entire spectral range 
(0.85-2.45 /xm). An identical instrumental configuration 
was used for the observations of field dwarfs of known 



spectral type (Testi et al. 2001 ; 2002b ). Data reduction 
and calibration was performed as described in Testi et 
al. ( p00l[ p002al ). 



2.3. Broad-band photometry 

On August 1 and 3 2001, we obtained moderately deep 
Gumi-i integrations using DFOSC and the Danish 1.54 m 
telescope at the ESO La Silla observatory. Following stan- 
dard bias, flat fielding and sky subtraction, typically 6 in- 
dividual 15 min dithered frames were coadded to produce 
the final images. Photometric calibration was a chieved by 
observing a set of stars from the Landolt (|l992|) catalogue, 
for which i-band AB magnitudes were computed using the 
transformations given in Fukugita et al. ( |1996| ). Given the 
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Table 1. Sample objects and i-band photometry 



(1) 

Object 
(ISO#) 




(2) (3) 
Coordinates 
(J2000.0) 




(4) 
i-band 
(AB mag) 


(5) 
Other 
Names 


023 


16 


26 


18.8 


-24 


26 


09 


20.34 ±0.15 


SKS1-10 


030 


16 


26 


21.4 


-24 


25 


59 


16.50 ±0.10 


SKS1-13 


















GY5 


032 


16 


26 


21.7 


-24 


44 


43 


16.26 ±0.10 




033 


16 


26 


22.2 


-24 


24 


05 


21.73 ±0.20 


SKS3-13 


















GY11 


102 


16 


27 


06.5 


-24 


41 


50 


15.75 ±0.10 


GY204 


160 


16 


27 


37.4 


-24 


17 


58 






164 


16 


27 


38.6 


-24 


38 


39 


18.18 ±0.10 


SKS1-49 


















GY310 


176 


16 


27 


46.3 


-24 


31 


41 




GY350 


193 


16 


28 


12.2 


-24 


11 


37 







References for C olumn 5. SKS: Strom et al. (1995); GY: Greene 
& Young (Ejgf) 



Table 2. Derived Object Properties 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Object 


ST 


T e ff 


L* 


A v 


M* 


(ISO#) 




(K) 


(Lo) 


(mag) 


(Uj) 


023 


M7 


2650 


0.04 


8.0 


30-50 


030 


M6 


2700 


0.07 


3.0 


40-80 


032 


M7.5 


2600 


0.06 


2.0 


30-50 


033 


M8.5 


2400 


0.008 


7.0 


8-12 


102 


M6 


2700 


0.08 


3.0 


40-80 


160 


M6 


2700 


0.04 


6.0 


30-60 


164 


M6 


2700 


0.09 


6.0 


40-80 


176 


M6 


2650 


0.07 


7.0 


30-70 


193 


M6 


2650 


0.1 


7.5 


40-80 



uncertainties in the transformations and the non perfect 
observing conditions, the uncertainties in the photometry 
are rather large. 

For all sources near-IR J,H,K S photometry is avail- 
able from the 2MASS second incremental data release. 
Additional L' and R-band photometry were taken from 
Comeron et al. ( |1998j ). 

3. Spectral classification and stellar properties 

The observed near-infrared spectra obtained at the TNG, 
normalized to the mean flux in the interval 1.1 - 1.75 /xm, 
are shown in Fig. |2| and || We derive for each object ef- 
fective temperature and luminosity in the following man- 
ner. We first obtain the extinction and spectral type by 
comparing the source spectra to those of field dwarfs. We 
then use the derived extinction value to obtain the effec- 
tive temperature through the comparison with reddened 



model atmospheres. The luminosity is computed from the 
dereddened J-band magnitude using the appropriate bolo- 
metric correction derived from the model atmosphere. 

The first step is illustrated in Figure ^, where we com- 
pare the observed spectra with a set of reddened field 
dwarfs, also obtained at the TNG with the same instru- 
mentation (Testi et al. |2001| , |2002b| ). We adopt the ex- 



tinction law appropriate for p Oph (R=4.2; Cardelli et 
al. |1989D . 

The overall shape of the spectrum from 1 to 2.4/im 
depends strongly on the spectral type of the object and 
extinction along the line of sight. There is, however, a 
degree of degeneracy, so that a cooler, less reddened object 
looks similar to a hotter, more reddened one. Therefore, 
we have also considered other features, such as the shape 
of the H band, the drop due to water absorption at the 
red edge of the J band and the intensity of some of the 
features visible in the J band. 

Although some of these characteristics depend some- 
what on the gravity, so that we cannot expect a perfect 
match between the young BDs and the field dwarfs, the fits 
are very good for most objects. Of the nine targets, three 
have extinction Ay< 3 mag, and six >6.0 mag; eight out 
of nine objects have a spectral type M6-M7, with ±one 
subclass uncertainty, while #033 has a later spectral type 
(M8.5). We give the results in Table 2, Columns 2 and 5. 
Note that even if the extinction is determined in the wave- 
length range 0.8-2.4 /zm, for convenience we express it in 
terms of Ay, the extinction in the visual, (Ay= Aj/0.313; 
Cardelli et al. |l989| ). 

In Fig. [| the same TNG spectra are compared to low 
gra vity, l ogg = 3.5, model stellar atmospheres (Allard et 
al. 2000| ; |200l| ), smoothed to the appropriate resolution 
and reddened using the value of Ay derived above. We 
obtain from this comparison the best value of the effec- 
tive temperature T e s, as well as a check on the adopted 
value of Ay. Our estimates of T e g have an uncertainty 
of typically ±100K; we assign values of 2600-2700 K to 
all the objects of spectral type M6-M7.5, while #033 is 
definitely cooler (T e g~ 2400). The vales of T e g are given 
in Table 2, Column 3. The robustness of these results 
and the uncertainties in the derivation of Ay and T e ff are 
described in more detail in Appendix [b| 

Since all the objects in our sample have a mid-infrared 
excess, we have considered the possibility that excess emis- 
sion is present also in the near-infrared, and is affecting our 
determination of the stellar parameters. We have there- 
fore subtracted the maximum contribution expected from 
an irradiated disk (flared, seen face-on; see §4) from the 
observed spectra. We found no significant change in the 
derived stellar parameters. 

The luminosity of the objects is shown in Table 2, 
Column 4. It has been computed from the dereddened 
J flux, and the ratio of the J to the total flux given by 
the appropriate stellar atmosphere model. These bolomet- 
ric cor rectio ns are virtually iden tical t o those of Wilking 
et al. dHH ) and Leggett et al. ([2002^ . For all stars, the 
adopted distance is D=150 pc. The uncertainties on L* 
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are difficult to determine accurately. We estimate that 
they probably range from 20% to 30%, mostly due to 
uncertainties on Ay- The bolometric correction for the 
J band changes very little with the atmospheric parame- 
ters, but an additional uncertainty (the same for all stars) 
may come from the uncertainty in the assumed distance. 

Finally, we have performed a last check on the reliabil- 
ity of our estimated parameters using our i-band photom- 
etry. For each star, we computed from model atmospheres 
and extinctions theoretical values of the magnitude in the 
i, J, H, K and L bands by convolving the flux distribu- 
tion with the appropriate filter responses. The results are 
shown in Fig.|], where we compare model predictions and 
observations for the 6 stars for which i-band photometry 
was obtained. [] The agreement of the i-band observed and 
predicted magnitudes is generally rather good, given the 
extreme sensitivity of the model predictions to the exact 
shape of the i-band filter, with the possible exception of 
#032, which would need Ay=3 mag, rather than the 2 
mag determined from the comparison with field dwarfs. 
The corresponding change in luminosity would be of 35%. 

Fig. U shows the location of the nine p Oph objects in 
the HR diagram. In the three panels, we overlay them to 
three different sets of e volut ionary tracks, co mpute d by 
D'Antona & Ma zzitel li (|l997|) , Chabrier et al. ( ^000| ) and 
Burrows et al. ( 199?|) , respectively. The derived masses 
(Table 2, Column 6) depend on the adopted tracks, hence 
we report the corresponding range of values. All objects 
appear to be very young, with ages lower than 1 Myr and 
probably of the order of a few 10 5 yr. It is well known 
that at such ages evolutionary tracks are not very reli- 
able (Baraffe et al. 2002 ), and that the parameters derived 
from the location on the HR diagram are only indicative. 
However, in spite of the uncertainties in both tracks and 
observations, we estimate that our sample contains one 
very low mass object (#033), with a mass of only ~8- 
12 Mj (Paper II), and a group of objects with masses in 
the BD range, of which about half (#023, #032, #160 
and #176) are very likely BDs. 

The clustering of eight out of nine of our objects in 
a narrow region of the HR diagram is a result of our se- 
lection criteria and can be understood as follows. TTS in 
p Oph have typical ages of 1 Myr, with very few stars as 
old as 3 Myr (Palla & Stahler |2000| ). The lack of older 
BDs in the sample is easy to understand, since the lim- 
ited sensitivity of the ISOCAM survey (especially at 14.3 
/im) strongly biases towards the highest luminosity, hence 
the youngest sources, and we expect to find in our sam- 
ple only BDs younger than the average TTS. Older, more 
massive objects could, in principle, fall in our sample. In 
practice, we found that this was not the case, given p Oph 
typical age. We have applied the procedure adopted by 
Bontemps et al. ( 2001 ) to a "theoretical" star with mass 
0.2 M Q and age of 2 Myr, using model-predicted J, H, K 



magnitudes (Baraffe et al. 1998) and Ay<9 mag; such a 
star would have a computed luminosity higher than our 
selected upper limit (L+< 0.04 L Q ), and would therefore 
not be included in our sample. Younger stars would be 
even more luminous. 

4. Disk models 

All nine objects have mid-infrared fluxes measured with 
ISOCAM in at least two b ands (centered at 6.7 and 14.3 
/im; Bontemps et al. 2001 ). In three cases, there are ad- 
ditional ISOCAM observations in three narrower bands, 
centered at 3.6, 4.5 and 6.0 /xm; Comeron et al. 
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1 Note that in Paper II, Fig. 4 shows in the inset broad-band 
fluxes of #033 (GY 11) dereddened by Ay=7.5 mag, rather 
than 7.0 as quoted. 



The ISOCAM points are shown for each object in Fig. 
together with our calibrated and de-reddened TNG spec- 
tra. 

For each system we compute the SED predicted by 
disk models, assuming that the disk is heated by the ra- 
diation of the central object. We ignore in this paper any 
possible viscous heating within the disk (see §6 for a brief 
discussion) . We follow as in Paper I and II the method out- 
lined by Chiang & Goldreich ( |l997j ; CG97), with s ome im- 
provements and modifications (Natta et al. |2001 ; Chiang 
et al. 2001 ). CG97 consider a disk in hydrostatic equilib- 
rium in the vertical direction (flared), and describe at each 
radius the vertical temperature structure of the disk in 
terms of two components: the disk surface, i.e., the exter- 
nal layer of the disk which is optically thin to the stellar 
radiation, and the disk midplane. These models allow a 
quick and reasonably accurate description of the expected 
SED, more than adequate for the purposes of this paper. 

The disk is a scaled-down version of TTS typical disks. 
It extends inward to the stellar radius, and outward to 
R D = 1 x 10 15 cm (67 AU). The total mass is M D ~0.03 
M*, and the surface density varies as R . The dust in the 
disk midplane h as op acity n = 0.01(A/1.3mm) _1 cm 2 g _1 
(Beckwith et al. 1990|) . For the dust on the disk surface, we 
take the mixture of carbonaceous materials and silicates 
that provides a good fit to the S EDs of several pre-main- 
sequence stars (Natta et al. 2001 ), i.e., a MRN distribution 
of graphite and astronomical silicates with dn/da oc a -3 5 , 

drain — 

ioo A, a ma x = l A* m ! 30% of cosmic C and all Si into 

grains. 

The results of the model calculations are shown in 
Fig. |(| The stellar parameters (T c ff, L*, M*) are taken 
from Table 2. As pointed out in Paper I, most of the 
disk parameters are irrelevant for the calculation of the 
mid-infrared disk emission, or appear in combinations, 
and cannot be determined individually (see also Chiang 
et al. 2001). As long as the disk midplane remains opti- 
cally thick to mid-infrared radiation, the only parameters 
that affect the SED in the near and mid-infrared are the 
geometrical shape of the disk (i.e., the flaring angle), the 
inclination to the line of sight and, to some degree, the 
disk inner radius Ri. There is also some dependence of 
the shape of the SED on the surface dust model; how- 
ever, since the luminosity intercepted and re-radiated by 
the optically thin surface layers is fixed, variations due 



Natta et al.: BD disks 



5 



- p-Oph ISO#023 


i i 1 i i i i _ 
A v =8.0 mag — 

-. 


" Ml 5 - 


_ 






M6.0 ✓ 




: / / : 




" M7.5 /'"* 




■ ^ A/ 




M8.5 







H 


1 1.5 


3 



4 - p-Oph IS0#030 

M4 5 




E 2 




Wavelength (jurn) 



1 1.5 2 

Wavelength (^m) 



__i i i I i i i i L_ 

1 1.5 2 

Wavelength (^m) 




1 1.5 2 

Wavelength (M m ) 



1 1.5 2 

Wavelength (M m ) 



E 2 




ft'-"" 



i i I i i , i L_ 

1.5 2 
Wavelength (/Ltm) 



4 - p-Oph IS0#16 



E 2 




Fig. 2. Observed TNG/Amici spectra of the sample objects. In each panel, we show (solid line) the spectrum of one 
obj ect com pared with the reddened spectra of field M-dwarfs ( dotted lines) of different spectral types, from Testi et 
al. ( 2002b ). All spectra are normalized to the mean flux in the 1.1-1.75 /xm range and shifted with constant offsets for 
clarity. The field dwarf spectra have been reddened by the value of Ay shown in each panel. 



to (reasonable) changes of the grain properties are well 
within the uncertainty of the existing observations. 

The upper solid curves in Fig. | show the SEDs of 
flared disks with i?,=R^, seen face-on. They all have 
strong silicate emission at 10 (im and a rather flat spectral 
slope between the two ISO bands at 6.7 and 14.3 /mi, of or- 
der a ~ 0.6 — 0.8 (vFy (xv a ). If, rather than extending all 
the way to the stellar surface the disk is truncated further 
out, as predicted by magnetospheric accretion models in 
TTS, at each radius the surface of a flared disk intercepts 



and reprocesses a larger fraction of the stellar radiation. 
The disk emission increases correspondingly at all wave- 
lengths but in the near-infrared, where one is sensitive to 
the lack of the hottest disk dust. A model with R4 ~ 3R* 
is shown (dashed line) for #033, where, as discussed in 
Paper II, the inner hole may account for the large ob- 
served mid-infrared excess. 

Large variations of the predicted SED occur if the disk 
shape changes. On each panel, we show the predictions 
of geometrically thin, "flat" disks (lower solid lines), i.e., 
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Fig. 3. Sam e as Fig. g, but in this case the red dotted spectra are reddened theoretical atmospheric models (Allard 
et al. 2000), with T e g as labelled and log g=3.5. 



disks where the grains are not well mixed with the gas, 
but have collapsed onto the disk midplane. Also for these 
models, we have adopted the CG97 formalism, which re- 
mains adequate in all the cases where the disk heating is 
dominated by the stellar irradiation. If the surface contri- 
bution to the SED were negligible, one would recover for 
these disks the well known temperature profile T oc i?~ 3 / 4 
and the power -law slope of the SED vF v oc v^l z (Adams 
and Shu 1986 ). Our calculations show that also in flat 



disks the surface contributes to the mid-infrared flux, as 
shown by the presence in the SED of the silicate feature 
in emission; however, the midplane emission is larger than 



the surface contribution at all wavelengths but in the re- 
gion ~ 8-12 ^m, where the silicate feature dominates, so 
that the spectral slope between the two ISO points is al- 
ways very close to 4/3. At all wavelengths larger than ~ 
2.5 /mi, the emission of a flat disk is significantly lower 
than that of a flared one. 

Finally, we show on three Panels of Fig. || the pre- 
dictions of tilted flared disks, seen by the observer with 
inclinations of 69° (#102), 80° (#164) and 86° (#193) 
respectively (0° for face-on disks). 

The comparison of the ISO observations to the model 
predictions shows that irradiated disk models can account 




Fig. 4. Comparison between dereddened observed magnitudes (i, J,H,K, L'), shown by squares, and the prediction of 
model atmospheres with parameters as in Table 2 (dashed lines) and of the same model atmospheres with additional 
disk emission (flared, face-on; solid lines) for the six stars for which i-band magnitudes are available. Disk models are 
described in §4. 



for the observed mid- infrared excess. More precisely, and 
in spite of the large uncertainties of the ISO data, inspec- 
tion of Fig. H shows that there are five stars out of nine 
(#030, #032, #102, #160, #176) that are extremely well 
fit by flat disk models. Two objects (#023 and #033) 
seem to require flared, face-on disks, while two others have 
a lower mid-infrared excess, consistent with disks seen 
rather edge-on. However, given the large error bars and 
the model uncertainties, most objects with flat disks are 
also consistent with flared disk models with large inclina- 
tion, as shown for the case of #102. 

5. Discussion 

5.1. Photospheric Parameters 

The agreement between the object spectra and those of 
field dwarfs is good beyond our expectations. The largest 
differences are of 20% at most, generally at the peak of 
the H band, with no systematic difference between ob- 
jects with large or low extinction, nor between stronger 
and weaker sources. Even the relatively narrow features 
that appear in the spectra around 1.1 /im are often well 



matched in the two sets of spectra. This indicates that, 
at the resolution of our observations, one should not ex- 
pect strong gravity effects. We have checked that this is 
indeed the case by comparing model atmosphere spectra 
smoothed to the obs erved resolution for stars of different 
gravity (Allard et al. |200l[ ), ranging from 3.5 to 6.0. All the 
models with gravity in the interval 3.5-5.0 are practically 
identical, at our spectral resolution and in this tempera- 
ture range. 

The comparison of our spectra with model atmosphere 
predictions is somewhat less satisfactory, especially in the 
H band, where the shape of the feature peaked at about 1.7 
/im (resulting from water absorption features at shorter 
and longer wavelengths) is narrower in the models than 
observed, and around 1.3 /im, where the models tend to 
predict more emission than is observed. Note, however, 
that this is not always the case (see, for example, #032 
and #164). Still, the agreement is in general rather good, 
with differences that never reach more 30%, again with no 
dependence on the extinction nor on the observed signal. 

The comparison of our determinations of the photo- 
spheric parameters of individual objects with previous 




Fig . 5. H R diagram for thre e sets of evolutionary tracks: D'Antona & Mazzitelli ( 1997|) in t he left panel, Chabrier et 
al. ( |2000| ) and Baraffe et al. ( |1998| ; for the 0.2 M ) in the mid panel, Burrows et al. ( |l997| ) in the right panel. Solid 
lines refer to objects of different mass, as labelled: hydrogen burning stars in black (red), deuterium burning BDs in 
medium grey (green), objects below the deuterium burning limit in light grey (cyan). Isochrones are shown as dotted 
lines, and labelled with the appropriate age. On each panel, the location of the nine observed objects is shown by dots 
with error bars. 



spectroscopic determinations in the literature shows that 
in some cases there is good agreement, while in others 
there are discrepancies that are not easily understood. 



For example, Wilking et al. (1999) assign similar spectral 
types to #023, #030, but a significantly later one (M8.5) 
to #164, based on K band R^300 spectroscopy. For the 
same object, Luhman & Rieke ( 1999] ) estimate a spec- 
tral type M7, based on intermediate resolution K band 
spectroscopy, similar to our classification M6. The same 
authors, on the other hand, attribute to #030 a somewhat 
earlier spectral type (M5-M6). The case of #033 (GY 11) 
has been discussed in detail in Paper II. A likely reason for 
differences in the spectral classification is that our scheme 
is based on the overall spectral shape, while the others 
rely on fitting individual spectral features, which in the in- 
frared s how la rge scatter for late M objects (e.g., Luhman 
& Rieke |1998|). 



On the more general issue of the effective temperature 
scale of young BDs, we attribute temperatures in the range 



2600-2700 K to our group of objects with spectral types 
M6-M7.5. Our only object with later spectral type (M8.5) 
has T e ff=2400±100 K. In a preliminary analysis of our 



sample field dwarfs (Tcsti et al. 2002b), we derive a similar 
effective temperature-spectral type correspondence. This 
is not significantly different from the scale used by Wilking 



et al. (1999) in their study of candidate BDs in p Oph. It 
is, however, at odds with some recent results, that tend to 
attribute to young BDs of similar spectral types temper- 
atures highe r than our values (Lucas et al. 2001; Lodieu 



2002J) . Further work, on larger samples of BDs in 



et al. 

young star forming regions is clearly required. 



5.2. The disk hypothesis 

The comparison between models and observations, dis- 
cussed in the previous section, proves that the mid- 
infrared excess associated to many young BDs can be ac- 
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Fig. 6. Disk and photosphere predicted SEDs. In each panel, the red dots with error bars show the ISOCAM observed 
fluxes (Comeron et al. 1998; Bontemps et al. 2001). The black solid line is the dereddened and calibrated TNG/Amici 
spectrum. The green jagged line shows the SED of the photosphere. The combined SED of the photosphere plus disk 
is shown by blue lines; in each panel, the two solid curves refer to face-on flared (upper curve) and flat disks (lower 
curve), with Ri=R*. For #033, the dot-dashed curve shows the SED of a face-on, flared disk with Ri=3R*. Finally, 
we show on three panels the SEDs of tilted flared disks (dotted lines), seen by the observer with inclinations of 69° 
(#102), 80° (#164) and 86° (#193) respectively (0° for face-on disks). 



counted for by the emission of circumstcllar disks heated 
by the radiation of the central object. 

Few disk properties are constrained by the existing ob- 
servations, and we do not want to overinterpret our re- 
sults, given the large uncertainties of the observed fluxes, 



and the simplicity of the adopted models. However, in 
our limited sample of nine stars we find disks of different 
flavours, and, in particular, an indication that many BDs 
may have flat disks. If we consider also the three objects 
in Cha I studied in Paper I, we have three objects with 
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clear evidence of flared disks, and nine where flat disks 
seem more appropriate, although we cannot rule out al- 
most edge-on flared disks for some of them (see also Apai 
et al. 2002| ). This is potentially an interesting result, since 
it seems natural to associate flat disks with dust sedi- 
mentation toward the midplane. In our selection of ISO 
sources, we have an obvious strong bias against objects 
with flat disks, since we required that the sources were de- 
tected by ISO in both bands. So, the fact that our objects 
with the lowest 6.7 fim fluxes (Cha Hal and #033) have 
flared disks is not surprising. However, there is no bias 
against selecting flared disk objects of higher luminosity, 
and we find only one (#023). The possibility of dust set- 
tling in these very young low-mass objects is intriguing. 
However, it needs to be confirmed by high-quality photo- 
metric observations at longer wavelengths, before entering 
into further speculations. 

The ejected embryos hypothesis does not exclude that 
BDs may have a small, and therefore short-lived, circum- 
stellar disk. Estimates by Bate et al. ( 2002| ) give disk radii 
of about 20 AU or less. The existing infrared data do 
not allow us to rule out such possibility, since the SED 
of a model with R D =20 AU will differ from the SED of 
a disk with Rd=75 AU only at wavelengths > 40 ^m. 
The mass of the disk is not predicted by the Bate et 
al. ( 2002| ) calculations, nor constrained by the existing ob- 
servations, since the only constraint we can set is that the 
disk has to be optically thick in the mid-infrared. This, 
however, only requires a disk mass of 10 -5 -10 -6 M©(or 
MD/Mstar^ 10~ 4 ), which is still consistent with a typical 
disk (having M D /Mstar~ 0.03, R D =75 AU), truncated at 
Rd=20 AU. Until far-infrared and millimeter data become 
available, the only way to validate these models is to de- 
termine the fraction of disks in unbiased samples of BDs 
of known age. 

Finally, one should remember that our analysis relies 
on the assumption that the ISO sources coincide with the 
objects we identify in the near-infrared. In some cases, 
this is likely to be true (see Appendix |a| and the discus- 
sion of #033 in Paper II). In other cases, it is impossible 
to check the validity of this assumption, given the large 
ISO beam and the presence of other red objects in the 
near-infrared images. However, the good agreement be- 
tween the observations and the model predictions, which 
depend essentially only on the stellar properties we derive 
from the spectroscopy, is encouraging. Further tests of the 
association of the observed mid-infrared excess with the 
identified stars could be obtained by accurate images in 
the L and M bands, where we predict that the disk emis- 
sion should be dominant (see Fig. ^|). 

6. Conclusions 

We have discussed in this paper a sample of nine very low- 
mass objects in the p Oph star forming region that have 
evidence for circumstcllar warm dust. We selected from 



and the 14.3 /xm bands, relatively low extinction and low 
luminosity. We determined first if these BD candidates 
were indeed bona-fide BDs, and then we checked if the ob- 
served infrared excess was consistent with the predictions 
of disk models, similar in properties to those associated to 
T Tauri stars. 

Our strategy was very successfull. The low-resolution 
near-infrared spectra obtained at the TNG allowed us to 
determine for each object spectral type and extinction, 
by comparison with field dwarfs observed with the same 
instrumental set-up, as well as effective temperature and 
luminosity, by comparison with model atmosphere predic- 
tions. The comparison with various sets of evolutionary 
tracks on the HR diagram shows that all the nine sources 
are very young, low-mass objects. In particular, one (#033 
or GY 11, already discussed in Paper II) has a mass of 8-12 
Mj, while the others have masses in the BD mass range; 
four of them are very likely bona-fide BDs. 

In all objects, the mid-infrared excess is consistent 
with the predictions of disks irradiated by the central ob- 
ject. We find no evidence of strong accretion occurring in 
these systems, based on the fact the observed near-infrared 
fluxes are dominated by the emission of the photospheres, 
and there is very little contribution (if any) from hot dust. 
However, it is not clear to which degree the near-infrared 
excess in very low-luminosity objects is a sensitive indi- 
cator of accretion (see, for an example of an actively ac- 
creting object with no near-infrared excess, Fernandez & 
Comeron 2001), and this issue should be explored more 



the ISOCAM sample of Bontemps et al. (|2001j) those ob- 
jects that have mid-infrared detections in both the 6.7 



quantitatively in the future. 

The existing data indicate that the disks must be op- 
tically thick at mid-infrared wavelengths; in some cases 
they must be flared (i.e., gas- rich with well- mixed dust 
and gas), while in others it is possible that they are geo- 
metrically flat, i.e., that dust has settled to the disk mid- 
plane. However, data at longer wavelengths are necessary 
to further investigate this point, and we do not want to 
put too much weight on this rather weak evidence. 

In the same sobering vein, we want to point out that 
our results do not discriminate yet between different for- 
mation mechanisms, namely between the possibility that 
BDs form from the gravitational collapse of individual, 
very low-mass cores, and the ejected embryo theory. We 
fit the observed mid-infrared excess with a scaled-down 
version of disks around the more massive TTS. This, how- 
ever, just implies that "normal" disks can account for the 
existing observations, since few parameters are actually 
constrained. As already pointed out in Paper I and II, 
only observations at long wavelengths can measure the 
disk radius and mass, since the lower limits that we can 
derive from the conditions that the disk is optically thick 
in the mid-infrared are hardly significant. 

Having stressed all the limitations of our results, let us 
now point out that this is the first sample of very low mass 
objects in a star forming regions where evidence for cir- 
cumstellar disks has been found and investigated in detail. 
Our accurate near-infrared spectroscopy, which allows us 
to estimate a reliable value of the mass of the objects, 
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proves that disks exist around low mass objects, well into 
the range of brown dwarfs. In one case, #033, our data 
provide strong indications that an object with mass close 
to or below the deuterium burning limit also has a cir- 
cumstellar disk. In addition to providing the beginning of 
a census of disk properties around BDs, our models indi- 
cate that the excess due to the cold disks irradiated by 
a central BD can only be detected by deep photometry 
in the L and M bands. We expect that major progress in 
our understanding of BD formation will be obtained by 
combining near-infrared low resolution spectroscopy with 
photometry in J,H,K,L,M of unbiased (i.e., not a-priori 
selected because they have a mid-infrared excess, as here) 
samples of BD candidates in star forming regions of dif- 
ferent age. 
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Appendix A: K.„ or K' finding charts 



Fig. A.l shows 2.2 /im finding charts for eight of the nine 



sources discussed in this paper. The grey scale images have 
been obtained at the TNG with NICS (#023, #030, #032) 
by us on July 2001, at ESO/UT1 with ISAAC (#033) 
and at ESO/NTT with SOFI (#102, #160, #164, #176). 
The ISAAC and SOFI data have been extracted from 
the ESO Science Archive; they have been originally ob- 
tained for ESO proposals 63.1-0691, 65.1-0576, 67.C-0325 
and 67.C-0349. Finding cha rt for #193 can be obtained 
from the 2MASS database at http:/ /irsa. ipac.caltech.edu. 
For sources #023, #030 and #033, contours show the 



ISOCAM-LW1 emission (Comeron et al. 1998). 



Appendix B: Accuracy in the determination of the 
photospheric parameters: Ay, 
Spectral Type, and T e g 

In this appendix we discuss the accuracy of the method 
used to derive extinction, spectral types and effective tem- 
peratures from the low-resolution near infrared spectra. 
For this purpose we will discuss the derivation of the pa- 
rameters for three extreme cases, the source with the high- 
est extinction (#023), the source with the lowest extinc- 
tion (#032), and the later spectral type source (#033). 

The procedure starts with the determination of the ex- 
tinction and spectral type based on the comparison with 
field dwarfs spectra obtained with the same instrumental 
setup (Testi et al. 2001 ; 2002b ). For each source we com- 
pare the observed spectrum with those of the field dwarfs 



reddened by different amounts, the value of the extinc- 
tion is varied until the best match is found with some of 
the field dwarfs. In practice, we varied the value of Ay 
in steps of 0.5 mags, since smaller variations cannot be 
significantly distinguished in the comparison. 

In the top panels of Figure B.l we show this compari- 
son for source #023 and for three values of Ay =7, 8, and 
9. The best matches, based on the shapes of the J and H 
bands, and the fact that the source spectrum at K cannot 
be lower than the field dwarfs, to allow for a possible con- 
tribution from disk emission, are found for Ay =8 mags 
and spectral types between M6 and M7.5. 

Once the spectral type and extinction have been esti- 
mated, the source spectrum is compared with appropriate 
surfac e grav ity (Log(g)=3.5) model atmospheres (Allard 
et al. 2000 ), which provide an estimate for the effective 
temperature and an additional che ck on the extinction. 
In the bottom panels of Figure B.l we present this com- 
parison for source #023. The comparison with the atmo- 
spheres confirms that the estimate of Ay =8 mags pro- 
vides the best match, and show that the two models that 
match more closely the observed spectrum are those with 
temperature 2600 and 2700 K, hence the best estimate 
that we derive is T e ff=2650 K, with an uncertainty of 
- 100 K. 



In Figure B.2 we show the procedure for the source in 
our sample with the lowest value of the extinction (#032). 
Given the lower extinction, the spectrum clearly shows 
some features in the J-band that correspond to the com- 
bination of several blended photospheric absorption lines 
and bands. Also in this case we show the comparison with 
field dwarfs for three values of Ay spaced by 1 magni- 
tude. In this case, it is possible to obtain a good match 
with a later spectral type and a lower extinction (Ay = 1 , 
M8.5); however, the match with M7.5 and Ay =2 mags is 
the best. This choice of the extinction is confirmed by the 
comparison w ith the model atmospheres in the bottom 
panels of Fig. B.2. The effective temperature is estimated 
to be 2600±100 K, as for the previous source. 

The source for which the derivation of the parame- 
ters is most complicated is #033 (Fig. BJ5) because of 
the high extinction and the lower signal to noise of the 
spectrum. In this case good matches can be obtained for 
the lower extinction and later spectral type, Ay =6 mag 
M9.0, as well as higher extinction and earlier spectral 
type, Ay =8 mag M7.5. The comparison with model at- 
mospheres would favor the higher extinction. Thus, also 
considering the uncertainties in the model atmospheres at 
these spectral types (see e.g. Leggett et al. 2001), our best 
estimates are: Ay =7.5±1 mag, M8.5 with one subclass 
uncertainty and T eff =2400±100 K. 
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Fig. A.l. Finding charts for eight of the nine sources: 2.2 /im images are shown as grey scale, ISOCAM-LW1 images 
are shown as contour plots. See text for details. 
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Fig. B.2. Same as Figure B.l, but for source #032, the values of the extinction in this case are: Ay =1, 2, and 3 mags. 



